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ABSTRACT 

Gravitational wave noise associated with unresolved binary stars in the Galaxy is studied 
with the special aim of determining the upper frequency at which it stops to contribute at 
the rms noise level of the proposed space-born interferometer (LISA). The upper limit to this 
background is derived from the statistics of SN la explosions, part of which can be triggered 
by binary white dwarf coalescences. The upper limiting frequency at which binary stochastic 
noise crosses LISA rms sensitivity is found to lie within the range ss 0.03 — 0.07 Hz, depending 
on the galactic binary white dwarf coalescence rate. To be reliably detectable by LISA, the 
energy density of relic cosmological background per logarithmic frequency interval should be 
tt GW h 2 W0 > 10~ 8 at / > 0.03 Hz. 

Subject headings: Binaries: close - gravitation - waves - white dwarfs 
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1. Introduction 

Binary systems constitute at least half stellar galactic population (<~ 10 11 stars) and are reliable sources 
of gravitational waves (GW). For a binary system consisting of two solar-mass stars, the characteristic 
time of orbital decay due to angular momentum removal by GW becomes shorther than the Hubble time 
(w 15 x 10 9 years) if its orbital period is less than about 14 hours. Observations of binary pulsars with 
a neutron star (NS) as secondary component (Taylor 1992) provides us with the strongest observational 
evidence for this fundamental process, and the most compact NS or black holes (BH) binaries merging 
due to GW emission are considered as primary real targets for the initial GW iterferometric LIGO-type 
detectors (Abramovici et al. 1992). An enormous energy is released during the coalescence of a compact 
NS+NS binary (typically, of the order of 10 53 ergs), but the galactic merging rate of NS+NS binaries is 
fairly small, 1CT 4 - a few 10~ 5 per year (see Lipunov, Postnov & Prokhorov 1997 for more detail), and to 
detect an acceptable number of such mergers per year (the experimentalists usually quote the number of 3 
events per year with the signal-to-noise ratio 3-5), they are planned to be studied from distances up to 200 
Mpc (e.g., Thorne 1987). Typical frequencies of GW emitted during compact binary NS coalescence are 
about 100-1000 Hz, depending on the mass of the stars involved. 

Orbital frequencies at which a binary system may be observable are limited by the size of the 
components; the evolution of the orbital separation strongly changes when one of the stars fills its critical 
Roche lobe. Typically, mass exchange between the components has much larger effect on the binary 
separation than the orbital angular momentum removal due to GW, so the orbital frequency of a given 
binary system may fall within a wide range from a fraction of hour to several years, depending on the initial 
parameters and details of evolution. The evolution of any star, however, ends up with the formation of a 
compact remnant (BH, NH or white dwarf (WD)), so when two such compact stars remains in a binary, its 
orbital evolution is totally controlled by the removal of orbital angular momentum by GW (the possible 
rare exceptions are very hot young WD with a strong stellar wind or the cases where a low-mass WD 
companion is evaporated by a strong pulsar emission; we will not consider them here). When only GW 
emission drives the orbital evolution, a simple analytical treatment is in order. 

As the number of binary stars in the Galaxy is very large, the GW emitted by them (at strictly twice 
the orbital frequency if the orbit eccentricity is zero) form a stochastic background in the frequency range 
from 10~ 7 to ~ 1 Hz (Mironovskij, 1965; Rosi and Zimmerman 1976; Lipunov and Postnov 1987; Lipunov 
ct al. 1987; Hills et al. 1990; Lipunov et al. 1995). Being interesting by itself, this background, however, is 
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viewed as a noise burying a possible cosmological gravitational wave background (CGWB), which bears the 
unique imprint of physical processes occurring at the very early (near-Plankian) age of the Universe (see 
e.g. Grishchuk 1988 for a review). 

A stochastic GW background is commonly measured in terms of the energy density per logarithmic 
frequency interval related to the critical energy density to close the Universe, £Igw = dEaw / d\n f / p cr c 2 
(p cr ~ 1.9 x 10~ 29 /i 2 00 g cm~ 3 where ft-ioo = Hq/WO km s _1 Mpc -1 is the present value of the Hubble 
constant, c is the speed of light). Cosmic microwave background fluctuations experimentally detected in the 
last years (Smoot 1992; Strukov et al. 1993) put upper bounds on CGWB of order Slew — 1CP 14 at both 
LISA and LIGO frequencies assuming equation of state p — — e at the inflationary stage of the Universe. 
In this model, CGWB spectrum is inversely proportional to the frequency keeping Q.gw{I) = constant at 
/ > 10~ 14 Hz (Rubakov et al. 1982). If real, this CGWB has no chance to be detected with ongoing GW 
detectors (see Schutz 1995, 1996; Allen 1996). We should note, however, that the measurements of cosmic 
microwave background temperature fluctuations put constraints on the CGWB only at very low frequencies 
Jh ~ Ho ~ 10~ 18 Hz, so the dependence of CGWB energy density upon frequency is very crucial. 

The situation, however, starts to change after a recent processing of the COBE data has provided us 
with some fresh information about the power-law spectral index of primordial perturbations (Bennett et 
al. 1996; Brukhanov et al. 1996). The point is that the spectrum of the primordial perturbations may 
be directly recalculated into the spectrum of relic GW background (see Grishchuk 1996 and references 
therein). According to Grishchuk (1996), the fact that the COBE data point to some deviations from 
Harrison-Zel'dovich spectrum of the initial perturbations translates into the deviations from p = — e 
equation of state at the epoch of inflation which would lead to serious changes in the CGWB as a whole. 
Now its energy density may increase with frequency reaching O-gw = 10~ 8 at / = 1CP 2 Hz. This conclusion 
holds regardless of the real cause of the observed CMB temperatrure fluctuations (cither they are due 
to density fluctuations or relic gravitational waves). If Grishchuk's calculations are correct, this opens a 
novel possibility of detecting CGWB with the proposed LISA space interferometer (Schutz 1996). The 
conventional technique for detection of the stochastic GW backgrounds assumes the cross-corelation of 
outputs of at least two independent interferometers (Grishchuk 1976; Compton & Schutz 1996 and reference 
therein). The LISA project does not plan to have two independent interferometers, however the strength of 
the signal (Qgw = 10~ 8 ) is predicted so high that probably the signal could be directly detected. 

This important finding raises the question: At which frequency GW confusion limit from binary 



- 5 - 



systems becomes lower than the detectability threshold dictated by the signal to noise ratio equal one? 
In other terms, beginning at what frequency can we be sure that no known noise sources of astrophysical 
origin exist and therefore if detected, only cosmological background can contribute at these frequencies? 

The aim of the present paper is to answer this question using new calculations of binary stellar 
evolution (see Lipunov et al. 1995, 1996b, 1997 for more detail and full references). The structure of the 
paper is as follows. In Section 2 we derive the upper limit on the binary stochastic background in our 
Galaxy at frequencies 1 mHz - 0.1 Hz (LISA diapason). At these frequencies, the background is mostly due 
to binary WD mergers. We use the SN la statistics to constrain galactic binary WD merging rate 1Z < 1/300 
yr _1 . The upper limit on the background at these frequencies is hu m (f ) ~ 4 x 10 _20 (//10 _3 Hz)~ 2 / 3 for 
the inverse-average distance to typical binary WD 10 kpc. In Section 3, we calculate the binary galactic 
background using Monte-Carlo modeling of binary star evolution in the Galaxy. The comparison of the 
calculated backgrounds with the proposed LISA sensitivity is given in Section 4. In Appendix A we give an 
alternative derivation of equation (Q) for the part of the background formed by coalescing binary stars only. 
In Appendix B we briefly discuss cosmological effects. 



2. Upper bounds on the binary stochastic GW background 

It is widely recognized that merging compact binary stars (WD+WD, NS+NS ...) determine the 
high-frequency part of binary stochastic background (e.g., Lipunov, Postnov & Prokhorov 1987; Hils, 
Bender & Webbink 1990). In this Section we wish to show that this portion is mostly significant in the 
LISA frequency band and provides an upper limit on this background in general. 



2.1. GW backgrounds formed by merging binary stars 

Consider a model galaxy consisting entirely of binary systems. Let us assume a stationary star 
formation rate (which describes well the situation in our Galaxy). We will be interested in only LISA 
frequency range, 10~ 4 — 10 _1 Hz, inside which only coalescing binary white dwarfs and binary neutron 
stars contribute. Even if binary neutron stars coalesce at a rate of 1/10000 yr in the Galaxy (Lipunov et al. 
1987; 1995; Tutukov and Yungelson 1993), their number still should be much smaller than the white dwarf 
binaries, and in this section we restrict ourselves to considering only binary WD. 

As mentioned in the Introduction, the stochastic GWB may be fully characterized by the energy 
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emitted per logarithmic frequency interval. In the case of a binary system this energy is exactly equal to 
the change in the orbital energy of two stars during the time period needed for the GW frequency (which 
is twice the orbital frequency for circular orbits) to pass the frequency interval A/ ss /. Therefore, in the 
stationary situation the net energy emitted in GW will be determined by the rate 1Z at which the systems 
enter the specific frequency interval and the rate of their orbital frequency change /. Clearly, if the binary 
orbit evolves only due to GW emission, the result will depend on the rate 1Z only. 

Indeed, the stationarity implies that the number of binary WD per unit logarithmic frequency interval 
may be determined from the continuity equation 

dN/d\nf = N(f)=Kx(f/f). (1) 

The total energy emitted per second per unit logarithmic frequency interval at / by all such binaries in the 
galaxy is 

dE/(dtd\nf) ee L(f) = £>;(/) = Hf)N(f) = L(f)K x (///) , (2) 

i 

where L(f) oc (Aif) 10 ^ 3 is the characterictic GW luminosity at frequency / which is dependent on the 
so-called "chirp mass" M. of the binary system: 

M = M(fi/M) 3 / 5 (3) 

(M is the total and fj, reduced masses). We use M. for some average mass of the typical binary (see 
discussion followed eq. below). For the orbital frequency change due to GW we find 

(f/fhw = (2/3)(E orb /E orb ) GW (4) 

(we have used the fact that E or {, oc M\Mij a and the third Kepler's law for binary semimajor axis a), so if 
E or b = {dE/dt)Q\Y we arrive at 

L(f) = (2/3)TZE orb (M,f). (5) 

For an isotropic background we have 

tt G w(f)PcrC 2 - L(f)/(4irc(r) 2 ) - 1lE orb (M, f)/(^c{r) 2 ) , (6) 

where (r) is the inverse-square average distance to the typical source. Strictly speaking, this distance (as 
well as the chirp mass) may be a function of frequency since the binaries characterized by different chirp 
masses M may be differently distributed in the galaxy. We are highly ignorant about the real distribution 
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of binaries in the galaxy, but taking the mean photometric distance for a spheroidal distribution in the form 

dN oc exp[— r/r ] exp[— (z/zq) 2 ] 

(r is the radial distance to the galactic center and z the hight above the galactic plane) with ro = 5 kpc and 
zq = 4.2 kpc with (r) w 7.89 kpc is sufficient for our purposes. 

Substituting E orb ~ Mc 2 (Mf) 2 ^ 3 into equation (j|) we obtain 

*WCf) « 2 x lO- 8 ^ 1O o(//lO- 3 Hz) 2 / 3 (7W/Af ) 5 / 3 ((r)/lOkpc)- 2 / l r o 2 o , (7) 

where TZmo — 72-/(0.01 yr _1 ) is the galactic rate of binary WD mergers. 

In terms of the characteristic dimensionless amplitude of the noise background that determines the 
signal-to-noise ratio when cross-correlating outputs of two independent interferometers (cf. Thorne 1987, 
eq. [65]) we have 

h c (f) = (1/27T) (ff //)C 

m 7.5 x lO- 2O 7e^o(//lO _3 Hz)- 2 / 3 (M/M ) 5 / 6 ((r)/lOkpc)- 1 

irrespective of Hq (naturally). Explicit derivation of equation (^) is given in the Appendix A. Equation (^) 
shows that at high frequencies of interest here the GW background is fully determined by the galactic rate 
of binary WD mergers and is independent of (complicated) details of binary evolution at lower frequencies 
(the examples of calculated spectra at lower frequencies see in Lipunov & Postnov 1987; Lipunov, Postnov 
k Prokhorov 1987; Hils et al. 1990, and below). 

All above considerations are valid for frequencies at which more than one binary system fall within the 
logarithmic frequency interval Af rs /. Formally this limiting frequency is specified by the requirement 
N(f) > 1, i.e. 

/ < f cr = 0.1Hz (M/M Q )- 5 / s TZ% s , (9) 

which is, however, already close to the limiting orbital frequency for a double WD binary. As the number of 
systems per logarithmic frequency interval rapidly decreases with frequency, N(f) oc /~ 8 ^ 3 , the stochastic 
background starts forming shortly below this limiting frequency. Note that combining equations (||) and 
(|) we can find the relationship h(f cr ) > 8 x 10" 22 (/ cr /0.01Hz) 2 / 3 (Al/Af o ) 5 ' /3 (r/10kpc)" 1 above which 
binary-merger-formed GW background appears; this limit depends on the mean galactic distance r, which 
is known much more reliably than binary WD merger rate 1Z. For M. < 1M Q , however, this boundary lies 
below the rms LISA sensitivity at all frequencies. 
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2.2. Notes on the extragalactic binary GW background 

The derivation of the binary background given above (see also Appendix A) becomes more accurate 
for extragalactic binaries. The contribution from extragalactic binaries may be shown to be smaller than 
from galactic ones at all frequencies (see, e.g. Hils et al. 1990, Lipunov et al. 1995), but for completeness 
we give here the final result in terms of Cl G w which bears a clear physical meaning. 

At the detector's frequency / = /'/(I + z )-, galaxies lying inside the proper volume dV(z) at redshift z 
contribute 

p cr c 2 dn GW (f) = L(f')/(47rcD(z) 2 )n(z)dV(z) (10) 

where L(f) is the GW luminosity at the source frequency / within each galaxy, D{z) is the luminosity 
distance, ng is the present density of galaxies and n(z) — n G (l + z) 3 is the density of galaxies at the 
redshift z without evolutionary effects. In the Universe with zero curvature dV{z) = 4ird 2 (z)d(d(z)) 1 where 
d(z) = D(z)/(1 + z) is the proper motion distance. Taking into account that the present density of galaxies 
can be rewritten through the fraction of luminous matter baryons, 

n G ps 0.013(Mp C - 3 )(f2 b /0.005)^ 00 ■ (H) 
equation (|l0|) can be recasted in the form 

dCl GW {f) = (2/3)Q b K {E orb (M, /(l + z))/M G c 2 ) (1 + z)d(d(z))/c (12) 

where M G is the baryon mass of a typical galaxy per which the merging rate 1Z is calculated. 



For flat (CI = 1) standard cosmological model without A-term d(z) = (2c/Hq)(1 — 1/y/T+z). Noticing 
that E or b(f(l + z)) scales as (1 + z) 2 / 3 (eq. [Al] from Appendix A), we may integrate equation ([l2] ) up to 



z* - the redshift at which first WD mergers had started - to obtain 

n G w(< z*) = {6/7)n b t H K(E orb /M G c 2 )((l + z») 7 / 6 - 1) 

= (S/7)Q b t H n(M/M G )(x/2) 2 ^((l + z,) 7 / 6 - 1) 

(13) 

fa 10~ 9 (fV0.005)ftioo 

(lO 11 X/M G )(X/Af ) 2 / 3 (//lHz) 2 / 3 / lu , 1 o((l + ^) 7/6 - 1) • 
This expression bears very clear physical meaning simply as the fraction of energy emitted in GW over the 
Hubble time (t H = 2/3H ) by WD mergers with respect to the rest- mass energy of baryons in stars in the 
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Universe. The ratio between the cosmological and galactic backgrounds is 

n GW {z*)/n GW = (r) 2 f \l + z) 2 ^n(z)dV{z)/D(z) 2 

J (14) 

w 0.05(fi b /0.005)((r)/10kpc) 2 / lu) 1 ((l + z*) 7 / 6 - 1) 

regardless of the poorly known merger rates, which shows that cosmological binary background is generally 
a few times smaller than the galactic one. Strong source evolution with redshift would increase this ratio, 
but at present we cannot calculate reliably these effects. 



2.3. SN la rate as the upper limit of binary WD merger rate 

The galactic merger rate of close binary WD is unknown. One possible way to recover it is searching 
for close white dwarf binaries. A recent study (Marsh et al. 1995), revealed a larger fraction of such systems 
than had previously been thought. Still, the statistics of such binaries in the Galaxy remains very poor. 

If coalescing binary WD are associated with SN la explosions, as proposed by Iben & Tutukov (1984) 
and further investigated by many authors (for a recent review of SN la progenitors see Branch et al. 
1995), their coalescence rate can be constrained using much more representative SN la statistics. Branch 
et al. (1995) concluded that coalescing CO-CO binary WD remain the most plausible candidates mostly 
contributing to the SN la explosions. The galactic rate of SN la is estimated 4 x 10~ 3 per year (Tamman 
et al. 1994; van den Bergh and McClure 1994), which is close to the calculated rate of CO-CO coalescences 
(<~ (1 — 3) x 1CT 3 ). The coalescence rate for He-CO WD and He- He WD (other possible progenitors of SN 
la) falls ten times short of that for CO-CO WD (Branch et al. 1995). 

We note further that binary WD mergers as main progenitors of the SN la explosions start falling out 
of favour in the last years (see especially critical studies by Nomoto et al. 1996). Evolutionary, double 
WD are formed through the common envelope stage (Iben and Tutukov 1984, Webbink 1984). In 1990s, 
with the calculation of new opacities, the possibility appeared for a WD, accreting even at a very high rate 
from the secondary companion in a binary system, to avoid the common envelope formation because of 
efficient stellar wind (Hachisu et al. 1996). This mechanism decreases the double WD formation rate and 
explains why the observed space density of binary degenerate dwarfs is smaller than derived from Iben and 
Tutukov's scenario. Therefore, we may conclude that the observed SN la rate provides a secure upper limit 
to the double WD merger rate regardless of the evolutionary considerations. 
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3. A model of the galactic binary GW background 

In this section we discuss the galactic GW background calculated with the Scenario Machine code 
for binary population synthesis in our Galaxy (see Lipunov, Postnov & Prokhorov 1995, 1996b for full 
description of the code). Briefly, we use a model galaxy with the total stellar mass of 10 11 M Q and half 
mass in binaries. We assume a constant star formation rate of 1 Mq per year which is a reasonable 
approximation for our Galaxy. Assuming the Salpeter mass distribution function dN/dM oc (M/M Q )~ 2 - 35 , 
this star formation rate requires the minimal mass of the star to be M m i n = 0.1 M Q in order to produce the 
total stellar mass of the Galaxy during the Hubble time of 15 billion years. The initial mass ratios of the 
modelled binaries q — M2/M1 are assumed to be uniformly distributed from to 1 (see, however, discussion 
in Lipunov, Postnov & Prokhorov 1996a). The important evolutionary parameter - the efficiency of the 
common envelope stage acE ~ was fixed at 1. Another crucial evolutionary parameter, the distribution of 
the kick velocity imparted to a neutron star at birth, is not very important here as we are interested mostly 
in low-mass stars (with initial masses Mi < 10 M Q ) which evolve ultimately into WDs. 

The position of the Sun in the outskirts of the Galaxy makes the binary stochastic background 
anisotropic (see Lipunov et al. 1995). The exact calculations of this background thus would require the 
knowledege of the spatial distribution of binary stars in the Galaxy, which is uncertain, so we restrict 
ourselves to consider the average amplitude of the stochastic background using the mean photometric 
distance in the Galaxy 7.9 kpc as explained in the previous section. 

The calculated background is shown in Fig. 1 by the thick solid line. For comparison, we plotted 
Bender et al's data (filled quadrangles). 

EDITOR: PLACE FIGURE HERE. 



4. Discussion and conclusions 

The upper frequency at which the stochastic GW background formed by galactic merging binaries 
becomes smaller than the detector's noise limit h rms (f) may be derived from equation (|J): 

f Um « (0.02 Hz) x (h rms /W~ 20 yV 2 1ll 00 (M/M G ) 5 / 4 (r/W kpc)" 3 / 2 
The upper limit (H) is plotted in Fig. 1 for different rates of binary WD mergers TvLioo = 1, 1/3, 1/10, 1/30 
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assuming the chirp mass M. » 0.52Mq (as for two CO white dwarfs with equal masses M\ = M2 = 0.6M©). 
These lines intersect the proposed LISA rms sensitivity at 

/> Sum .« 0.03 - 0.07Hz . (15) 

This means that at frequencies higher than 0.07 Hz no continuous GW backgrounds of galactic origin 
are presently known to contribute above the rms-level of LISA space laser interferometer. The contribution 
from extragalactic binaries is still lower regardless of the poorly known binary WD merging rate (at least 
in the limit of no strong source evolution with z). Other possible sources could be extragalactic massive 
BH binary systems (e.g. Hils and Bender 1995). Their number in the Universe can be fairly high (e.g. Rees 
1997), but no reliable estimates of their contribution are available at present. The lower limit ( |l5| ) is already 
close to the LISA sensitivity limit at 0.1 Hz, but we stress that the assumptions used in its derivation are 
upper limits, so the actual frequency beyond which no binary stochastic backgrounds contribute may be 
three times lower. This precise limit depends on the details of binary WD formation and evolution which 
are still poorly known. 

Fig. 1 demonstrates that the calculated GW background intersects LISA sensitivity curve at frequencies 
~ 0.05 Hz, and Bender et al's curve at even lower frequencies ~ 0.01 Hz. The latter is probably due to 
Bender et al's curve being derived from observational estimate of double WD galactic density in the solar 
neighborhoods; we stress once more that once formed, the binary WD will evolve until the less massive 
companion fills its Roche lobe; unless the mass ratio is sufficiently far from one (cf. Webbink 1984), the 
merger should occur. Therefore Bender et al's curve provides a secure lower limit to the galactic binary 
stochastic GW background. We also note that if the coalescence of two WD with inequal masses is 
prevented by mass transfer process (then the orbital period of the system begins increasing), some features 
on the shape of the background may emerge (cf. discussion of the minimum orbital period for cataclysmic 
variables in Lipunov, Postnov & Prokhorov 1987). 

Presently, we cannot rule out the high galactic double WD merger rate (1/300-1/1000 yr ), and 
therefore can consider Slim to lie within the frequency range 0.01 — 0.07 Hz. We conclude that no GW 
background of galactic origin above this frequencies should contribute at the rms-noise level of LISA 
interferometer, and hence the detection of an isotropic stochastic signal at frequencies 0.03 — 0.1 Hz with 
an appreciable signal-to-noise level (which possibly may be done using one interferometer) would strongly 
indicate its cosmological origin. To be detectable by LISA, the power of relic GW background should be 
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^Gwhi Q > 10 8 in this frequency range. 

The work was partially supported by Russian Fund for Basic Research through Grant No 95-02-06053-a 
and by INTAS Grant No 93-3364. We thank Prof. Leonid Grishchuk for stimulating discussions. KAP 
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No 496057. 

A. Explicit derivation of equation (H|) 

Equation (@) can be derived explicitly, without using the notion of flaw (instead, one can use h c to 
determine VLqw] both means are, of course, fully equivalent). We reproduce here this equation keeping 
proportionalities 

h c (f) oc {l/r)ll 1/2 M 5/(i f- 2/3 . 
We find it convenient to introduce new variables: 

R g = 2GA4/c 2 , gravitational radius of the chirp mass , 
Ri = c/(jrf), light cylinder radius for orbital frequency , nf/ 2 
and their dimcnsionless ratio 

x(M,f)=x = (R g /R l )(xMf. 

This units allows us to write all relevant quantities in a compact form with clear physical meaning keeping 
G and c and facilitating numerical estimates. In these units, the orbital energy of a binary system with 
chirp mass Ai is 

E orb = -(l/2)GMiM 2 /a = -{l/2)Mc 2 {x/2) 2 l 3 . (Al) 



The energy flux per unit time per unit area carried by a GW is 

dE/{dAdt) = c 3 /(l6irG)(h + (t) 2 + h x (t) 2 ) 



(A2) 
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where h + and h x are wave polarizations, the overbar means averaging over several cycles of the wave. 
Switching to the frequency domain in usual way (see, e.g., Thorne 1987), we obtain 



dE/(r 2 dndf) = (c7G)(7T/V2)(|M/)| 2 + IM/)| 2 ) , (A3) 
where r is the distance to the source, dVl = dA/r 2 is the elementary solid angle. 

Let us describe the stochastic backgrounds formed by many independent sources within unit logarithmic 



frequency interval at / by the characteristic strain amplitude defined as (cf. eq. [A2]) 

N(f) 

hl(f) = 4G/(™ 3 )(1// 2 ) £ dEt/tfdndt) (A4) 

i 

where N(f) — TZf/f is the number of sources assuming the stationary source birth rate TZ (cf. eq. |]). We 
wish to show that this equation gives exactly the same result as equation (0). 



Substituting equations (A2) and ( |A3j ) into equation (A4), using equation (|lj), and introducing, as 
before, some effective distance and the chirp mass of a typical source (now we will not mark them by special 
symbols), we arrive at 

^(/)=^/<|M/)| 2 ), (A5) 

where (\h(f)\ 2 ) = (J2i=i l^«(/)| 2 ) stands for angle-averaged Fourier-components (with dimension [Hz -2 ]) of 
h(t): 

(\h(f)\ 2 ) = (1/4tt)(1// 2 ) / 2G/(7rc 3 )(dE/r 2 dfdn)dn 



(A6) 

= (l/47rr 2 )2G/(irc 3 )(dE/df)/f 2 . 
(Note that eq. ]A5[ can also be derived directly from Parseval's theorem and the expression for the variance 
per logarithmic frequency interval of a stochastic noise associated with unresolved independent sources). 

Now notice that for circular orbits (our initial assumption) dE/df = {dE or ^/ df) if far from coalescence 
(which is the case under consideration), so 

dE/df = -(l/3)(Mc 2 /f)(x/2) 2 ^ 3 ex M 5/3 r 1/3 , (A7) 



and from equation (A6) 

(|M/)| 2 ) = {l/l2^r 2 )R 2 /f 2 (x/2f' 3 (x M 5 ' 3 r V3 (A8) 
(cf. eq. [44] in Thorne 1987). Finally, from equation ( |A5| ) we obtain 

h 2 c (f) = (^/6 7 r)(i? 3 /r) 2 (x/2)- 1 / 3 // cx M 5/3 f~ i/3 , (A9) 



which yields exactly equation 
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Fig. 1. — Galactic binary GW background h c as given by Bender (1996) (filled quadrangles) and calculated 
for a model spiral galaxy with the total stellar mass 10 11 M Q (the solid curve). Average photometric distance 
7.9 kpc is assumed. Thin straight lines marked with 100, 300, 1000, 3000 are the analytical upper limit (eq. 
H) for binary WD merger rates 1/100, 1/300, 1/1000, and 1/3000 yr _1 in a model spiral galaxy, respectively, 
assuming M. = 0.52 M©. Straight dashed lines labeled by 10~ 8 , 10 -6 show GW backgrounds corresponding 
to constant Qgw- The proposed LISA rms noise level (h rms ) and sensitivity to bursts hsB — 5y5h rm s are 
also reproduced (cf. Thorne 1995; Fig. 14). 



